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We present the mechanism of the dynamical supersymmetry (SUSY) breaking at the 
metastable vacuum in the J\f = 1 U(N) SUSY gauge theory with adjoint superficlds. 
, The dynamical SUSY breaking is triggered by the non-vanishing D-term coupled to the 

Cl+ observable sector, and is realized by the self-consistent Hartree-Fock approximation of the 

NJL type while it eventually brings us the non- vanishing F-term as well. We numerically 
check the local stability of our metastable vacuum. 
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1. Introduction 
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Spontaneous breaking of SUSY occurs much less frequent compared with that of 
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internal symmetry in quantum field theory and has attracted much interest!^ of 
theorists for over the three decades. Mass hierarchy in elementary particle physics 
indicates that it is most desirable to break AT — 1 SUSY dynamically. In fact, the 
non-renormalization theore nP protects the generation of holomorphic operator in 
perturbation theory and instanton generated nonperturbative superpotentials have 
been the major source of dynamical SUSY breaking (DSB). 

We focus our attention on general TV" = 1 theory in four dimensions consist- 
ing of vector superfields and chiral superfields in the adjoint representation with a 
non-canonical gauge kinetic function. It has recently been shown in refsP that, in 
this general situation, SUSY is dynamically broken in the metastable vacuum. The 
mechanism that triggers the DSB is the condensate of the Dirac bilinear, forcing 
one of the order parameters D of SUSY to be non- vanishing. This is very much rem- 
iniscent of the Nambu-Jona Lasinio (NJL) theory^ of broken chiral symmetry and 
hence the BCS superconductivity^ being formulated in terms of the effective action 
of the auxiliary field whose stationary value is the order parameter. The method 
of approximation employed is the self-consistent Hartree-Fock (HF) approximation 
where the tree and the one-loop contributions are regarded as comparable. Once this 
mechanism operates, non- vanishing i^-term is shown to be induced and contributes, 
for instance, to the mass of the fermions. The mechanism requires massive adjoint 
scalars, in particular, the scalar gluons and, together with the feature that the D 
term triggers the breaking, is quite distinct from the previous proposals^' of DSB 
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both from theoretical and experimental perspectives. The overall U{1) where the 
non- vanishing D and the Nambu-Goldstone fermion (NGF) reside serves as the hid- 
den sector and no messenger field is necessary^ as non- vanishing third prepotential 
derivatives connect the U(l) sector with the observable SU(N) sectorP 

In the next section, we review the original reasoning that has led us to the D- 
term triggered DSB. We set up the background field formalism to be used in the 
subsequent sections, separating the three kinds of background from the fluctuations. 
In section three, we elaborate upon our treatment of the effective potential with the 
three kinds of background fields as well as the point of the HF approximation in refs. 
[3] In section four, we present a qualitative argument of our variational analyses of 
the effective potential. Treating F-term as an induced perturbation, we demonstrate 
that the stationary values (_D* , 93* , 93* ) are determined by the intersection of the two 
real curves. Numerical analysis is provided that demonstrates the existence of such 
solution as well as the self-consistency of our analysis. The second variation of the 
scalar potential is computed and the local stability of the vacuum is shown from 
the numerical data. Summary is given in the last section. 



2. The action, assumptions and some properties 

The lagrangian we consider is 

C = J d 4 9K{<S> a ,$ a ) + {gauging) + J d 2 9lm^T ab {^ a )W aa W b a 

+ (y d 2 8W{$ a ) +c.c)j . (1) 

At the lowest order in perturbation theory, there is no source which gives a vacuum 
expectation value (vev) to the auxiliary field D°: (Z?°)trcc = 0. The U{N) gaugino 
is massless at the tree level while the fermionic partner of the scalar gluon receives 
the tree level mass m a = m = (g OQ dodoW)tree- 

It is useful to summarize here a set of assumptions. 

(1) a general JV — 1 SUSY action of chiral supcrficld <& a in the adjoint representa- 
tion and the vector superfield V a with the Kahler potential K{& a , $ a ) with its 
gauging, the gauge kinetic superfield T Q {,(<I> a ) following from the second deriva- 
tives of a generic holomorphic function J r ($ a ), and the superpotential IU($ a ). 

(2) third derivatives of J r ($ a ) at the scalar vev's are non- vanishing. 

(3) the superpotential preserves N = 1 SUSY at tree level. 

(4) the gauge group is U{N) and the vacuum is taken to be in the unbroken phase 
oiU{N). 

In refs.[3l it was shown that the vacuum develops a non- vanishing vev of an auxiliary 
field D° in the HF approximation. The theory, therefore, realizes the D-term DSB. 
The relatively simple estimate has shown that the vacuum can be made long lived. 
Let us recall a few more key aspects. 
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The part of the lagrangian which produces the fermion mass matrix of size 2N 

is 

(-4^1^) (£)+<->■ (2) 

It was observed that the auxiliary D a field, which is an order parameter of JV = 
1 SUSY, couples to the fermionic (but not bosonic) bilinears through the third 
prepotential derivatives: the non-vanishing vev of D° immediately gives a Dirac 
mass to the fermions. Equation of motion for Z?-term implies 

(D°) = -^(9°° (J r 0cd^ d X c + JW d A c )}, (3) 

telling us that the condensation of the Dirac bilinear is responsible for (D°) 7^ 0. 
We diagonalize the holomorphic part of the mass matrix: 

MFa =\-^ aa D°) {dada W) J' ^ 

Note that the non-vanishing third prepotential derivatives are J-baa where a refers 
to the generators of the unbroken gauge group. The two eigenvalues of eq. ((!]) for 
each generator are 

Agi = (d a d a W)\£i Ag> = \ (l ± v/l + AL) , A^r = ^fjy - (5) 

It was also shown in refs. [3] that the non- vanishing F° term is induced by the 
consistency of our procedure of computation. (See also ref. [9]). This is because the 
stationary value of the scalar fields gets shifted upon the variation (the vacuum 
condition). The final mass formula for the SU(N) fermions is to be read off from 

4it 0) - -\(gva, a )(F a )rr + l -(T Qaa )(F n )\ a x a - l -{d a d a wwr 

+ ^(^aaW^(D°). (6) 

The mass matrix is read off as 

-^9 aa ToaaF , -^^/g™(ImT)™T 0aa D° 



M a = 



~^^g aa {lmT) aa FvaaD°, g aa d a d a W + g aa g a : al 

(7) 

We parametrize this matrix such that, in the case of F° = F° = 0, its form reduces 
to that of refs. [3J The quantities with multiple indices such as To aa receive U(N) 
invariant expectation values: (J-Q aa ) = (J'ooo) e.t.c. We suppress the indices as we 
work with the unbroken U(N) phase in this paper. The two eigenvalues of the 
holomorphic mass matrix are written as 



1 / /, / 1 



A<±> = (trX)A(±), A(±) = - l±W(l + z/) 2 + (l + -/) A 2 ) (8) 
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where A = —2m\^/m^, f = 2im\\/trM. 

3. The effective potential in the HF approximation 

In the HF approximation, one begins with considering the situation where one-loop 
corrections in the original expansion in h become large and are comparable to the 
tree contribution. In this section, we start the analysis of this kind for our effective 
potential. There are three constant background fields as arguments of the effective 
potential: <p = ip° (complex), U(N) invariant background scalar, D = D° (real) and 
F = F a (complex). The latter two are the order parameters of Af — 1 SUSY. 
Let us denote our effective potential by V. 

V = Utrcc + V c . t . + V x 

— loop • (9) 

The first term is the tree contributions, the second one is the SUSY counterterm 
and the last one is the one-loop contributions, which are explicitly given by 

I 



V tr 



-gFF - -(ImF")D 2 - FW' 



FW 



V c . t . = -ilm J d 2 9AW 0a W 0a = -i(ImA)L> 2 , 
A(e, 7 ) (|A(+)| 4 + |A(-)| 4 
-|A(+)| 4 log|A(+)| 2 -|A(-)| 4 log|A(-)| 2 + 



Ul-loop 



N 2 \trM\ 4 
32rr 2 



trM 



m s 


4 


m s 


4" 


trM 


log 


trM 





(10) 

(11) 

(12) 



where the scalar gluon mass is m s (ip, ip) = g 1 (ip, ip)W'(ip). The one-loop potential 
is calculated by the dimensional reduction scheme and the divergence reside in 



s 1 1 „ d 

A(e, 7 ) = -- 7+ -, e = 2--. 



(13) 



The divergence is subtracted by a counterterm associated with ImJ 7 " by setting up 
a renormalization condition 

1 d 2 V 



N 2 (dDf 



= 2c. 



(14) 



Note that this condition is set up at D = and the stationary point of the scalar. 



4. Stationary conditions, gap equation and stability of our false 
vacuum 



Now we turn to our variational problem. It is stated as 
dV 
dD 
dV 

of 

dV 

dip 



= 0, 

= and its complex conjugate, 
= and its complex conjugate. 



(15) 
(16) 
(17) 
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We will regard the solution to be obtained by considering eqs. (fT5j) and (IT71) first 
and solving D and ip for F and F: 

D = D*(F,F), ip = tp m (F,F), (p = (p*{F,F). (18) 

Eq. (HH) is then 

dV(D = D*(F,F),ip = ip*(F,F),(p = (p* (F, F),F,F 



dF 



= (19) 

-D,(/p,V?,-F fixed 



and its complex conjugate. These will determine F = F*, F = F*. 

Here, we are going to work in the region where the strength \F*\ is small and 
can be treated perturbatively. This means that, in the leading order, the problem 
posed by eq. (fPo)) and eq. (fTT|) becomes 

dV(D,<p,<p,F = 0, F = 0) =Q) 



dD 

dV(D, ip, <p, F = 0, F = 0) 9F(£», ^,^,^ = 0,^ = 0) 



= (21) 



and this problem does not involve the tree potential eq. (fTU|) except the D 2 term, 
as F and F are set zero. Eq. ([2TJ|) is nothing but the gap equation, while eq. (|2"Tj) is 
the stationary conditions for the scalar. This is the variational problem which we 
analyze. A set of stationary values (-D*, tp*,<p*) is determined as the solution. 

In this paper, we do not write down the explicit expressions for the gap equation 
and the stationary conditions for the adjoint scalar because of their complicated 
forms (See the last paper in refs. [3] for details) and the limitation of the space to 
describe. Instead, we give a qualitative argument how these conditions are solved 
and numerical results in the simplified case. 



Fig. 1. The schematic picture of the intersection of the two curves which represent the solution 
to the gap equation (the red one) and the tp flat condition (the blue one). The horizontal axis is 
denoted by tp/M and the vertical one by Ao- The values at the stationary point (Ao*,<£* = <£*) 
are read off from the intersection point. 



The stationary values (D* , tp* , tp*) are determined by the gap equation and the 
stationary condition. The solution to the stationary condition in the An profile is 
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determined as the point of intersection of the potential with the quadratic term 
having (p = (p dependent coefficients. Actually, it is a real curve in the full (Ao, ip — 
(p) plane. Likewise, the solution to the gap equation provides us with another real 
curve in the (Ao,iy5 = (p) plane. The values (Ao*,<y9* = ip*) are the intersection 
of these two. The schematic figure of the intersection is displayed in Figure [TJ By 
tuning our original input functions, it is possible to arrange such intersection. 

We study some numerical solutions to the gap equation and the stationary con- 
dition for ip in the real Ao case. In order to find ip* explicitly, the forms of the 
prepotential T and the superpotential W must be specified. Here, we take a simple 
prepotential and a superpotential as 

T= tr</? 2 + tros 3 , W= trip + — — trip 3 , (22) 

2N y 3\MN N 3\N V ' 

where m, M carry dimensions. In particular, M is a cutoff scale of the theory. 

Some numerical solutions to the gap equation and the stationary condition for 
tp in some parameter points are listed in the Table below. In these examples, we 
have taken some values of — f — . ... and to just for an illustration and the ratio 

Im(z+A) J 

|D,/D*| is evaluated. We can find that the D-term is smaller than the D-term in 
these examples. Also, the scalar gluon mass squared are positiv^f], which ensures 
a local stability of our vacuum. As a summary of our understanding, a schematic 



Ao* 




\F*/D*\ 


m% 


0.477 


0.707 (10000) 


0.524 (to < M) 


0.4998 


1.3623 
1.3623 
1.3623 


0.8639 (2000) 
0.5464 (5000) 
0.3863 (10000) 


0.224 (to < M) 
0.142 (to < M) 
0.100 (to < M) 


0.7463 
0.2986 
0.1492 



figure is drawn in Fig. [21 which illustrates the local stability of the scalar potential 
at the vacuum of dynamically broken SUSY in comparison with the well-known 
NJL potential. 

5. Summary 

We have proposed a new mechanism of DSB at the metastable vacuum in a M = 1 
SUSY U(N) gauge theory with an adjoint chiral super field and non-canonical gauge 
kinetic term, in which non- vanishing D-term vev is developed in the self-consistent 
HF approximation of NJL type. Non-zero D-term is also induced by such a D-term 
as well. We analyzed a gap equation for D-term and the stationary conditions for 
ip and D-term in the case that D-term can be treated as a perturbation comparing 



a In the Table, only the leading term of the scalar gluon mass is shown. See the last paper in refs. 
[3|for more details. 
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Vnjl 






D 



Fig. 2. Comparison of V sca i ar around the stationary value (D*,ip t ) with Vnjl- 

to _D-term. The solutions are obtained as the intersection point of two curves repre- 
senting the gap equation and the stationary conditions. Numerical examples of the 
solutions are found in the simplified case and the local stability of vacua we found 
is also confirmed by the numerical data. 
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